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Abstract. 

We report systematic measurements using 59 Co nuclear quadrupole resonance (NQR) 
technique on the cobalt oxide superconductors Naa;Co02T.3H20 over a wide Na 
content range x=0.25~ 0.34. We find that T c increases with decreasing x but reaches to 
a plateau for x <0.28. In the sample with x ^0.26, the spin-lattice relaxation rate 1/T\ 
shows a T 3 variation below T c and down to T ~ T c /6, which unambiguously indicates 
the presence of line nodes in the superconducting (SC) gap function. However, for 
larger or smaller x, 1/Ti deviates from the T 3 variation below T ~ 2 K even though 
the T c (^4.7 K) is similar, which suggests an unusual evolution of the SC state. In 
the normal state, the spin correlations at a finite wave vector become stronger upon 
decreasing x, and the density of states at the Fermi level increases with decreasing x, 
which can be understood in terms of a single-orbital picture suggested on the basis of 
LDA calculation. 



PACS numbers: 74.25.Nf, 74.70.-b 
1. Introduction 

The layered cobalt oxides Na x Co02 have attracted much attention in recent years, 
because of a rich variety of orders competing for the ground state. For 2=0.75, the 
material shows a large thermoelectric power [1,2], while the x=0.5 compound is a charge 
ordered insulator with Tco=50 K [3]. The magnetic properties are different for different 
Na contents [3-5]. At the side with x larger than 0.5, the DC susceptibility shows a 
Curie- Weiss temperature (T) variation, while at the low- x side the susceptibility is T- 
independent [3]. Most interestingly, superconductivity with T c ~4.5 K emerges around 
x=0.3 when water molecules are intercalated in between the Co layers [6]. It has been 
reported that T c can be tuned by varying x in the range of 0.25 < x <0.38 [7]. 
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One of the unsettled issues on the superconductivity is the symmetry of the gap 
function. Theoretically, a chiral(d + id) superconducting state was first proposed on 
the basis of RVB scenario [8-10]. In a previous study, we found that the spin-lattice 
relaxation rate 1/Ti shows no coherence peak below T c , and follows a T n variation with 
n— 2.2~1 at low T [11]. This suggests that the superconductivity is of non-BCS type, 
but it has been pointed out that the data are insufficient to distinguish between some 
exotic states [12-14]. It has been shown that both chiral and pure ci-wave states in the 
presence of disorder [12,13], or an odd-frequency s-wave state [14] can be compatible 
with such T- variation of T x . Since identification of the gap function is the first step 
toward the understanding of the superconductivity, more work is needed. 

Another issue is how the electron correlations in the normal state evolve with 
doping, which should help understanding the mechanism for the superconductivity. 
In the octahedral crystal electric field environment, two e g orbitals of Co are located 
much higher in energy above ti g orbitals. The local density approximation (LDA) band 
calculation shows that [15-17], among the three t2 g orbitals, the a\ g (3rf 32 2_ r 2, where 
^-axis is the threefold rotation axis) orbit lies higher in energy than the remaining two 
e g orbitals. If one takes this picture [8,9], then the Co 4+ state without Na has one 
electron in the a lg orbital (spin 1/2), with the e g orbitals fully occupied. This single- 
orbital situation resembles the undoped cuprates; adding Na dopes electrons into the 
ai g orbital. In contrast, multiple-band models were also proposed [18,19], which draw 
similarities to the manganese oxides, another important class of strongly correlated 
materials. 

In this letter, we address the above two issues by studying the evolution of the 
system with doping using NQR (nuclear quadrupole resonance) technique. We present 
intensive data on the superconducting and the normal states over a wide range of Na 
content. 

2. Experimental Results and Discussion 

The Na^CoCVyH^C) powder [20] was synthesized following Ref. [6]. The precursor of 
Nao^CoC^ was immersed in the Br 2 /CH 3 CN solution to de-intercalate Na. The amount 
of Br 2 is listed in Table 1, where X means the theoretical value necessary to remove 
all the Na. The Na content of the final product was determined by the ion coupled 
plasma method. X-ray diffraction indicates that the samples are single phase except 
the 5X (x ~0.34) and 10X (x ~0.32) samples in which minor secondary phases are 
present. The c-axis elongates as the Na content decreases, in agreement with previous 
report [7]. T c was determined from the ac susceptibility and the 1/Ti measurements. 
The previously reported sample [11] was made with 10X Br amount in a separate run 
and turned out to have an x ~0.31. All other samples were made in the same run. 
59 Co NQR measurements were carried out by using a phase-coherent spectrometer. The 
nuclear magnetization decay curve is excellently fitted to the theoretical curves [24], 
with a unique Ti component. 
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Table 1. Characterizations of the samples. 



Br amount 


Na content x 


c-axis length (A) 


a-axis length (A) 


Tc (K) 


5X 


0.34 


19.50 


2.822 


2.65 


10X 


0.32 


19.62 


2.822 


3.60 


10X [11] 


0.31 


19.59 


2.819 


3.70 


20X 


0.28 


19.73 


2.817 


4.75 


40X 


0.25 


19.80 


2.817 


4.70 


100X 


0.26 


19.79 


2.815 


4.60 



Figure 1(a) shows a typical example of the three NQR transition lines arising from 
the nuclear spin 1=7/2 of 59 Co. Figure 1(b) shows the ± 3/2«>±5/2 transition line 
for various Na contents. Firstly, the peak frequency increases as x decreases down to 
x=0.28. Analysis of three transition lines allows us to find that the vq increases as 
the Na content decreases down to x=0.28, but the asymmetry parameter r\ decreases 
monotonically with decreasing Na content, as seen in Fig. 1(c) and (d), respectively. 
Here vq and r] are defined as u Q = v z = 2 i(2i-i)h e ' 2 Qwz'i V = \ v x ~ v y \/v z , with 
Q being the nuclear quadrupole moment and §^r(ct = x,y,z) being the electric field 
gradient at the position of the nucleus [21]. The decrease of vq with increasing x is 
consistent with doping electrons into the a ig orbit [23]. Secondly, for the 5X (x ~0.34) 
and 10X (x ~0.32) samples, there appear multiple peaks for the transition. However, 
the additional peaks with smaller amplitude are not due to the minor impurity phase 
detected in the X-ray chart. Rather, we find that 1/Ti measured in different peaks 
differs only by 10% in magnitude and shows the identical temperature variation; a 
clear decreases was seen below T c . These peaks may have arisen from the ordering of 
Na ions [22] which are located above and below Co. More detailed discussion will be 
published elsewhere [23]. We note in passing that the multiple peaks seen in the ± 
3/2<->±5/2 transition are not resolvable in the ± 5/2^->±7/2 transition. 

Figure 2 shows the temperature dependence of the quantity \jT{T for various x. All 
measurements were done at the ± 3/2^±5/2 transition. For the 5X and 10X samples, 
1/Ti was measured at the main peak shown in Fig. 1(b). Two trends can be seen in Fig. 
2. First, the absolute value increases with decreasing x. Second, \jT{T increases upon 
lowering T, indicating the importance of spin correlations; the smaller the x the steeper 
the increase, indicating that the spin correlations become stronger as x is reduced. 

What is the nature of the spin correlations? One good way to tell is through 
inspection of the Knight shift, which measures the spin susceptibility and, unlike 
the bulk susceptibility, is not affected by the presence of magnetic impurities. In a 
superconducting single crystal with x ~0.3, we found that the shift is T-invariant from 
T c up to T=120 K [30], which is in good agreement with the report of Ning et al [31]. 
These results rule out the possibility of ferromagnetic spin fluctuations with wave vector 
q=Q, since in that case the Knight shift would follow the temperature dependence of 
1/TiT because both of them would be dominated by the susceptibility at wave vector 
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Figure 1. (a) 59 Co NQR spectra for the sample of x ~0.28. (b) The line shape of the 
±3/2 <-> ±5/2 transition for various Na contents, (c) and (d) Na content dependencies 
of the asymmetric parameter r\ and vq, respectively. 




Figure 2. The quantity \fT\T as a function of temperature for various Na contents. 
The curves are fits to 1/T X T = {l/T x T) + C/{T + 0), with resulting parameter C= 75, 
77, 111 and 123 Sec" 1 for a;=0.34, 0.32, 0.26 and 0.25, respectively. The values of 
(l/TiT) and 6 are plotted in Fig. 3. 



q ~ 0. 
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The results thus indicate that the spin fluctuations peak at finite wave vector. 
The most likely candidate is antiferromagnetic (AF) correlations [8-10,32], although its 
detailed g-dependence awaits clarification by future experiments, in particular neutron 
scattering. Below we analyze our data in terms of AF correlations, in an attempt to 
quantify the discussion. In a general form, 1/71 T is written as 1/71 T = J2 q AgX"(q, <jj)/uj 
in the uj — >-0 limit, where A q is the hyperfine coupling constant, and is the 

imaginary part of the g-dependent, dynamical susceptibility. If one assumes that there 
is a peak around the AF wave vector Q, then one may have the following approximation: 
1/71 T = (1/7\T) AF + (1/717%. Here (1/717% denotes the contribution from wave 
vectors other than Q, which usually is dominated by the uniform susceptibility xo- For 
the AF contribution, we use the formula for a two-dimensional nearly antiferromagnet, 
namely, (1/T\T)af — fq^ (Ref. [33]). We find that such model fits the data quite 
successfully, as depicted by the curves in Fig. 2. The resulting fitting parameters 
are shown in Figure 3. Note that a small value of 6 means that a system is close 
to the AF instability (0=0). The enhancement of the spin correlation (decrease of 
6) with decreasing x, as is the case in doped copper oxide high-T c superconductors, 
can be naturally understood in the single-orbital picture [8, 9] described earlier. In 
contrast, the multiple-band theory predicts Na content-insensitive, ferromagnetic spin 
correlations [19]. 

Turning to the quantity (1/71T) , its increase with decreasing x can be attributed 
to the increase of density of state (DOS) at the Fermi level. This result is consistent with 
the LDA band calculation which found that the dominant Fermi surface is a hole-like 
one and that the DOS increases as the electron number is reduced [16,34]. This aspect is 
different from the case of doped cuprates where the DOS decreases as decreasing carrier 
density. We summarize our results in Fig. 3. As x decreases from 0.34, T c increases, but 
reaches to a plateau of ~4.7 K for x <0.28. Our results for high value of x are in good 
agreement with the report by Schaak et al [7] . A plateau has also been reported recently 
by Chen et al and Milne et al [25,26]. It is interesting to note that as the parameter 
9 decreases and C increases, namely, as the electron correlation becomes stronger, T c 
increases. 

Next, we turn to the superconducting state. As can be seen already in Fig. 2, no 
coherence peak was seen just below T c in all samples, irrespective of Na content and/or 
T c value. However, the low-T behavior is sample dependent. We first show in Fig. 4 the 
results for the sample with x ~0.26. It can be seen that 1/71 follows a T 3 variation down 
to the lowest temperature measured (~ T c /6). This result indicates unambiguously that 
there exists line-nodes in the gap function. As we noted in previous publications [11,27], 
the gap function with line nodes generates an energy (J^)-linear DOS at low E which 
results in a T n (n—3) dependence of 1/71. in the previous sample, 1/71 did not follow 
an exact T 3 variation and became proportional to T below T ~ T c /2, [11]. It has 
been proposed that many exotic superconducting states can be compatible with the 
data [12-14]. Calculations incorporating impurity/disorder scattering showed that a 
d x 2_ y 2-\-id xy state which has a full gap is more consistent with the previous data [12]. 
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Figure 3. Various parameters obtained in this study plotted as a function of Na 
content, x. The curves are guides to the eyes. 



Our present result completely rules out the d x 2_ y 2+id xy state proposed by the t— J model 
in an isotropic triangular lattice [8-10,28]. Also, the odd frequency state proposed [14] 
is gapless [35,36], and is incompatible with our finding in the x=0.26 sample. However, 
if anisotropy of the t term is incorporated, the theories could be compatible with our 
experimental results. In fact, it has been recently shown that a small anisotropy in t 
(10~20%) lifts the degeneracy between d x 2_ y 2 and d xy states, leading to a stabilization 
of a pure d-wave state [29]. 

In other samples with either larger or smaller x, however, \jT\ deviates from the 
T 3 relation below T=2.0 K. To see this more clearly, we plot the T-dependence of 1/T\ 
for all samples in Fig. 5. In the inset of the figure is shown the normalized 1/Ti by its 
value at T = T c against the reduced temperature. What is responsible for the deviation 
from the T 3 variation in samples with x either larger or smaller than 0.26? Disorder 
or impurity is an unlikely candidate, since the samples with x ~0.28 and x ~0.25 have 
very similar or even higher T c and NQR line-width as the x ~ 0.26 sample. Clearly, 
this unprecedented property deserves further investigation, both experimentally and 
theoretically. Meanwhile, low-lying excitations due to nearby competing orders may be 
a possible candidate. For larger x, it is known that charge ordering [3,22] occurs. The Na 
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Figure 4. Temperature dependence of 1/Ti for the sample of x=0.26. The arrow 
indicates T c . The straight line depicts the T 3 variation. 



order could also give rise to the anisotropy of the hopping integral t mentioned earlier. 
For smaller x, on the other hand, it has been proposed that different charge/magnetic 
order may occur at £=1/4 [37]. 

3. Conclusion 

In conclusion, we have presented intensive 59 Co- NQR measurements on the cobalt 
oxide superconductors Na x Co0 2 -1.3H 2 over a wide Na content range x=0.25~ 0.34. 
We find that T c increases with decreasing x but reaches to a plateau for x <0.28. In the 
superconducting state, the spin-lattice relaxation rate 1/Ti shows a T 3 variation down 
to T ~ T c /6 for x ~0.26, which indicates unambiguously the presence of line nodes in 
the superconducting gap function. For both larger or smaller x, however, 1/Ti deviates 
from the T 3 dependence even through the T c value is similar. This unprecedented 
finding shows the richness of the physics of this new class of superconductors, and 
any pertinent theory must explain this peculiar property. In the normal state, the 
temperature dependence of \jT{T indicates that the spin correlations at finite wave- 
vector become stronger upon decreasing x, while the DOS increases. This aspect is 
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Figure 5. A comparison of the temperature dependence of 1/Ti in the 
superconducting state for various samples with different Na content. The inset shows 
the normalized 1/Ti by its value at T = T c against the reduced temperature. 



different from the case of doped cuprates but can be understood in terms of a single- 
orbital picture suggested on the basis of LDA calculation. 
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